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Abstract. Wc consider a system of diffusion processes that interact through their empirical 
mean and have a stabilizing force acting on each of them, corresponding to a bistable potential. 
There are three parameters that characterize the system: the strength of the intrinsic stabilization, 
the strength of the external random perturbations, and the degree of cooperation or interaction 
between them. The latter is the rate of mean reversion of each component to the empirical mean of 
the system. We interpret this model in the context of systemic risk and analyze in detail the effect of 
cooperation between the components, that is, the rate of mean reversion. We show that in a certain 
regime of parameters increasing cooperation tends to increase the stability of the individual agents 
but it also increases the overall or systemic risk. We use the theory of large deviations of diffusions 
interacting through their mean field. 
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1. Introduction. Systemic risk is the risk that in an interconnected system of 
agents that can fail individually, a large number of them fails simultaneously or nearly 
so, leading to the overall failure of the system. It is a property of the interconnected 
system as a whole, and not only of the individual components, in the sense that 
assessment of the risk of individual failure alone cannot provide an assessment of the 
systemic risk. The interconnectivity of the agents, its form and evolution, play an 
essential role in systemic risk assessment [BJ. 

In this paper we consider a simple model of interacting agents for which systemic 
risk can be assessed analytically in some interesting cases. Each agent can be in 
one of two states, a normal and a failed one, and it can undergo transitions between 
them. We assume that the dynamic evolution of each agent has the following features. 
First, there is an intrinsic stabilization mechanism that tends to keep the agents near 
the normal state. Second, there are external destabilizing forces that tend to push 
away from the normal state and are modeled by stochastic processes. Third, there 
is cooperation among the agents that acts as individual stabilizer by diversification. 
This means that in such a system we expect that there is a decrease in the risk of 
destabilization or " failure" for each agent because of the cooperation or diversification. 
What is less obvious is the effect of cooperation on the overall or system's risk, which 
can be defined in a precise way for the model considered here. We show in this 
paper that for the models under consideration and in a certain regime of parameters, 
the systemic risk increases with increasing cooperation. The aim of this paper is to 
analyze this tradeoff between individual risk and systemic risk for a class of interacting 
systems subject to failure. 

Perhaps a simple mathematical model of interacting agents having the features we 
want is a system of stochastic differential equations with mean-field interaction. Let 
Xj (t) be the state of risk of agent or component j, taking real values. For j = 1, . . . , N , 
the Xj (t) 's are modeled as continuous-time stochastic processes satisfying the system 
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of Ito stochastic differential equations: 

dxj{t) = -hU(xj(t))dt + 6(x{t) - Xj(t))dt + adwj(t), (1.1) 

with given initial conditions. Here —hU(y) = —hV'(y) is the restoring force, V is 
a potential which we assume has two stable states, and {wj{t)}^ =l are independent, 
standard Brownian motions. The parameter h controls the level of intrinsic stabi- 
lization and a is the strength of the destabilizing random forces. The interaction 
or cooperation is the mean reversion term with rate of mean reversion 9 and with 
x(t) := jf yij^-i x i{t) denoting the empirical mean of the processes, that is, the em- 
pirical mean of the individual risks. For 6 > the individual risk processes tend to 
mean-revert to their empirical mean, which is a simple but non-trivial form of co- 
operation. We take the empirical mean x(t) to be a measure of the systemic risk. 
The bi-stable-state structure of V(y) determines the normal and failed states of the 
agents. We will assume in this paper that U (y) — y 3 — y, so that V(y) — \y A — ^y 2 +c 
and we take c = since it is inessential. The two stable states are then ±1 and we 
let —1 be the normal state and +1 to be the failed state. The potential V{y) ensures 
that each risk variable Xj(t) stays around —1 (normal) or +1 (failed). The evolution 
of the system is characterized by the initial conditions, the three parameters (h, 6, a) 
and by the system size N. 

We have chosen a mean-field interaction because it is a simple form of cooperative 
behavior. More elaborate models are considered in SectionEl where some heterogene- 
ity is introduced between the components of the system. For mean-field models a 
natural measure of systemic risk is the transition probability of the empirical mean 
x(t) from the normal state to the failed state. More precisely, the mathematical prob- 
lem we address here is this: For N large we calculate approximately such transition 
probabilities and analyze how they depend on h, a and 9, the three parameters of the 
system. We are interested in a regime of these parameters for which there are two 
collective, that is, large N, equilibria centered around the normal and failed states. 
These two equilibria can be identified through the mean-field limit of the system, that 
is, the weak limit in probability of the empirical density of the agents risk Xj. Mean 
field models with multiple stable points, not only bistable ones, could be considered 
but their analysis is more involved while the main result about systemic risk, and 
dependence on the parameters (h, 0, a) and by the system size N, is clearly seen in 
the bistable model that we consider here. 

The mathematical analysis of bistable mean field models like (jl.ip was initiated 
by Dawson [51 118). including the mean field limit, the existence of multiple equilibria, 
and a fluctuation theory. Non-equilibrium statistical mechanics and phase transitions 
have been studied extensively in the sciences |19j . The large deviation theory that we 
use here was developed by Dawson and Gartner [T0l[TT] . In particular, they introduced 
and analyzed the rate function for large deviations associated with Ijl.ljl when N is 
large and with more general potentials |llj . Their theory may be considered as an 
infinite dimensional extension of the Freidlin-Wentzell theory of large deviations for 
stochastic differential equations with small noise [TBI HI]- The main result in this 
paper is the analysis of this rate function for small h. That is, for a shallow two- 
well potential, where transitions from one well (quasi-equilibrium) to the other are 
exponentially small in N, the "constant" in the exponent is small when h is small. 
Other mean field models have been studied in [33l [TBI EH El EEl [30l HI] , and large 
deviations results for various models can be found in [T21 [TJ ESI EH E21 HI H] • In [3 a 
general large deviations theory is developed for a model with both drift and volatility 
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interactions, as well as with degenerate noise, using weak convergence and optimal 
control methods. 

The main contribution of the paper as far as systemic risk theory is concerned is 
the demonstration that, within the range of the bistable mean field model (ll.lj) . while 
cooperation between agents decreases the individual risk of each agent, the systemic 
or overall risk is increased. This is discussed in detail in Section 16.41 in terms of the 
three parameters Qi,6,o~), with h small. The fact that reducing individual risk by 
cooperation or diversification can lead to increased systemic risk has been anticipated 
in macroeconomics and elsewhere and it has been extensively discussed, modeled, and 
analyzed in [31] [4] [SO] [TT1 [26l [32j [5] [3] [2TJ [24] . However, the dynamic phase transitions 
formulation and the large deviations theory exploited in this paper have not been used 
in the economics literature, to our knowledge. The use of coupled stochastic equations 
for modeling evolution of individual risk and the effects of interactions among agents 
is also considered in [5] 123) where there is some discussion regarding the economic 
interpretation of the variables {xj(t)}. They could, for example, represent some form 
of equity ratio in a very simple model in insurance or banking. 

The paper is organized as follows. In Section (2] we briefly review the classical 
mean- field limit in [5], and we discuss the intrinsic stability of equilibria [5] when h 
is small. Section [3] generalizes (jl.ip by replacing the rate of mean reversion 9 by an 
agent-dependent 9j . The mean- field limit and the explicit conditions are also studied. 
In Section 01 we carry out numerical simulations of both the homogeneous and the 
heterogeneous model in various parameter ranges. Section [5] uses the large deviation 
principle in [10] to formulate the dynamic phase transition of interest here, that is, 
the system transition from the normal state to the failed state. In Section [6j we 
specialize the large deviations theory when h is small so as to obtain a result from 
which the systemic risk as a function the basic parameters (h, a) can be assessed 
and interpreted. In Section [7] we introduce a formal expansion of the rate function for 
small h and obtain a reduced variational principle for the systemic risk that appears 
to come from a large deviations principle for an one-dimensional dynamical system. 
It gives, of course, the same results about systemic risk as described in Section [6] 
In Section [3] we discuss the case where there is diversity in mean reversion and it is 
shown that under some natural conditions the heterogeneous model is systemically 
more unstable than the homogeneous one. The technical details of the proofs are in 
the appendices. 

2. The Mean-Field Limit. We briefly review the mean field limit in [HI [TB] 
and carry out a small h analysis of results since they will be used in calculating large 
deviation probabilities. We want to analyze the systemic behavior of the interacting 
diffusion processes (jl.ll) . through their empirical mean x(t), but this is not possible 
in a direct way since is nonlinear. We consider instead the empirical density of 
Xj(t), which is a measure valued process that has a limit as N — > oo. Let Afi(R) be 
the space of probability measures endowed with the weak (Prohorov) topology and let 
C([0, T], Afi(R)) be the space of continuous Mi (R)-valued processes on [0, T] endowed 
with the corresponding weak topology. Define the empirical probability measure pro- 
cess X N (t, dy) := ^E^Li^ W^y) and note that Xjv € C([0, T], Mi(R)). The mean 
field limit theorem for X^, proved in [HI [EH], is as follows: 

Theorem 2.1. (Dawson, 1983) Assume that the force is U(y) = y 3 — y and 
that Xn(0) converges weakly to a probability measure vq. Then the measure valued 
process Xn converges weakly in law as N — > oo to a deterministic process with density 
u(t,y)dy € C([0,T], Mi(R)) ; which is the unique weak solution of the Fokker-Planck 
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equation: 

d 3 



with initial condition vq. 

By Theorem 12. 1[ we can analyze u and view Xjy as a perturbation of u for TV 
large. We may consider x{t) in the same way because x(t) = J yX]y(t,dy). However, 
the limit problem is infinitely dimensional, as is expected. 

Explicit solutions of (|2.ip are not available in general, but we can find equilibrium 
solutions. Assuming that £ = lim t _ ! . 00 J yu(t, y)dy, then an equilibrium solution u| 
satisfies 

^-f)«a-^[«-«)«|]+^ u ?= ' 



and has the form 



uf (y) = )== exp | - h^V(y) \ , (2.2) 



Z^2tt% I *5S 
with the normalization constant: 

w^t- 4 ^ -*>»}*■ 

Now £ must satisfy the compatibility or consistency condition: 

e = m(0 := j yu\(y)dy. (2.3) 

Finding equilibrium solutions has thus been reduced to finding solutions of this equa- 
tion. 

For U(y) — y 3, — y, £ = is a solution for (|2.3|) . With the same U(y), it can be 
shown (see also [H Theorem 3.3.1 and 3.3.2]) that there are two additional non-zero 
solutions ±£b if and only if ^m(O) > 1, and for given h and 8, there exists a critical 

a c (h,9) > such that ^m(O) > 1 if and only if a < a c (h,6). 

An explanation for this bifurcation at equilibrium is that when a > cr c , ran- 
domness dominates the interaction among the components, i.e., dixit) — Xj{t))dt is 
negligible. In this case, the system behaves like N independent diffusions and hence, 
by the symmetry of V{y), at any given time roughly half of them stay around — 1 
and half around +1 so the average is 0. When, however, a < <r c , then the interac- 
tive force is significantly larger (now o~dwj(t) is less important). Therefore all agents 
stay around the same place (either — & or +£&) and the zero average equilibrium is 
unstable. Since we want to model systemic risk phenomena, we assume that a < o~ c 
throughout this paper, and we regard — as the normal state of the system and 
+£& as the failed state. The calculation of transitions probabilities between these two 
states is our objective. 

For small h we can approximate the solution of (I2.3[) to order O(h) as follows. 

Proposition 2.2. For small h, the critical value a c can be expanded as 

<T c = <M + 0{h). (2.4) 
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In addition, the non-zero solutions ±£f, are 



±6 = ±\/i-3- 



1 - 2{a 2 /29) 
29 J l-3(a 2 /29) 



0(h 2 ). 



(2.5) 



Proof. See Appendix [A] □ 

From Proposition 12. 2[ we see the relation between the existence of the bi-stable 
states and the ratio a 2 /29: For a given 9, and for small h, (I2.3|) has non-zero solutions 
if and only if 3a 2 /29 < 1. Moreover, these non-zero solutions ±£t are generally not 
±1 since the magnitude |^| is less than 1. Note that the coefficient of order h in the 
expansion (|2.5|) depends significantly on and a. Thus, when 3a 2 /29 tends to 1, £j, 



in (12. 5p will not go to +00 while, in fact, ^ goes to 0. From the O(l) term in (|2.5j) . 
we also see that is roughly decreasing as a 2 /29 is increasing. 

3. Diversity of Sensitivities. We can generalize by allowing for agent 

dependent coefficients. We consider a particular case in which each agent can have a 
different rate of mean reversion to the empirical mean, that is, for j = 1, . . . , N, 



dxj 



V(xj)dt + adwj + 6j(x — Xj)dt, 



(3.1) 



and as before V{y) = 



— \y 2 . We consider the case where 8i,...,9 N take K 
distinct positive numbers, 9i, . . . , Qk- We define Ii = {j : 9j = 9;}, pi = \Ii\/N and 
■^N = pTF Sjeii Assuming that limAr^oo pi exists and is positive for all I, the 
limit of (Xh, . ■ . , X¥j) as N — > oo are the weak solutions (u\ , . . . , uk) of the set of K 
coupled Fokker-Planck equations. 

Theorem 3.1. Assume that U(y) — y a — y and that (Xjf(0), . . . ,X§(0)) con- 
verge weakly in probability to the probability measures (v 1 , . . . , v K ). Then the measure 
valued vector process (-Xjy, ■ ■ ■ , Xjy) converges weakly as N —> oo to the weak solution 
(til, . . . , uk) of the system of the Fokker-Planck eguations: 



~ot 



dy- 



■ui - 9i 



d_ 

dy 



K 



y^2piui(t,y)dy ■ 



1=1 



ui 



d 

+ h—[U(y) Ul } 
ay 



(3.2) 



dt 



1 d 2 

UK = -7T^ U K - 9 A '- 



2 dy 2 



i^2piui(t,y)dy 



1=1 



UK 



h—[U{y)u K ] 



with initial condition [y , . . . , V s ). 

Proof. See Appendix IB.ll for the outline of the proof following [18] . □ 
The equilibrium solutions {wfjj/Ei have the form 



ut ti (y) 



1 



: exp ■ 



«V ""267 

1 



2^ o 2 



(3.3) 



J 26, 



: exp 



2tt 



\j£-h-%V(y))d y , 



20, 



• 26, 



() 



and £ must satisfy the compatibility condition 



£ = MO -=^2 PI / yut s (y)dy 



(3.4) 



1=1 



For U(y) = y 3 — y, £ = is the trivial solution of (|3.4[) . and a simple extension of 
Theorem 3.3.1 in |9], shows that there are two sets of non-trivial solutions {uf 
and {u'f _£ i }iL 1 if and only if ^m(O) > 1. The numerical simulations presented in 
the next section show that diversity in the rate of mean reversion can have significant 
impact on the stability of the mean-field model. 

As in the homogeneous case, we can get an approximate condition for equilibrium 
bifurcation for small h. 

Proposition 3.2. The compatibility condition {3.$ has non-zero solutions if 
and only if a < crf w . For small h, erf™ has the expansion 



Proof. See Appendix lB~2l □ 

We note that diversity does affect the threshold condition and makes the analysis 
more difficult. The non-zero solutions ±£b can be computed approximately when h is 
small: 



Higher order terms in the expansion of (13.51) can also be obtained but we will omit 
them in this paper. In the following Proposition we show that a^ lv < a^ om °, where 
^homo _ a ^ £ ne cr }ti ca [ value (|2.4I) of the homogeneous case. 

Proposition 3.3. With 6 = Y,f =1 Pl&l, we have a h c omo > af v for small h. 
Proof. See Appendix IB.3I □ 

This result shows that when there is diversity the parameter region of existence 
of equilibria ±£b is smaller than in the homogeneous case . From this observation we 
can anticipate that these equilibria are less stable in the presence of diversity, and 
this is confirmed next by numerical simulations and analytically. 

By noting that ^ omo = ^/l - (a 2 /a c homo ) 2 + 0(h) and £ h div = y/l - [a 2 /af v ) 2 + 
0(h), we have the following corollary: 

Corollary 3.4. With 9 = J2?=i Pl®l, we have 1 > $ om0 > for small h. 

4. Numerical Simulations. Before going into a detailed analysis of the models, 
we carry out numerical simulations of and (|3.1[) so as to get a quick impression 
of their behavior. We discretize with a uniform time grid, and let X™ denote the 
simulated Xj at time nAt. 

4.1. Homogeneous Model. We simulate (|1.1| using the Euler scheme 





(3.5) 




n+l 



A]' - hU(X™)&t + aAW? +1 + 0(— Y, X k~ Xf)At. 



(4.1) 
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N=100, h=0.1, o=1, 6-6 N=100, h=0.1, <J=2, 6=6 N=1Q0, h=0.1, o=4, 6=6 

dm(0)/d^=1 .01 26, 3tj 2 /26=0.25 dm(0)/d^=1 .0002, 3c 2 /26=1 dm(0)/d^=0.957D2, 3tj 2 /20=4 




Fig. 4.1. Simulations for different cr. The system has two (statistically) stable equilibria when 
a is below the critical value or otherwise has single stable state 0. For small h, 3<r 2 /20 < 1 is the 
approximate criterion. 



N=1Q0, [1=0.1,0=1,9=5 N=100, h=0.1, o=1, 0=1.5 N=100, h=0.1, c=1, 6=0.5 

dm(0)/d^=1.0141, 3o 2 /26=0.3 dm(0)/d^=1.0024, 3c 2 /26=1 dm(0)/d^=0.81904, 3o 2 /26=3 




Fig. 4.2. Simulations for different 8. The system has two stable equilibria if 6 is above the 
critical value or otherwise has single stable state 0. For small h, 3o 2 /29 < 1 is the approximate 
criterion. 



We take U(y) =y 3 -y, = l 1 X° = -1, At = 0.02, and let {AWf}^ n be independent 
Gaussian random variables with mean zero and variance At. In the figures presented, 
the dashed lines show the numerical solutions of the compatibility equation (|2.3[) . 
£ = As noted earlier, if ^rn(O) < 1, then = m(0) is the unique solution and 

is a stable state. Therefore we should observe that the systemic risk fluctuates around 
0. If ^m(O) > 1, there are two additional non-zero solutions ±£(, = m(±£ b ) and ±£ b 
are stable while is unstable. We also know that when h is small, the condition 
^■m(O) > 1 can be simplified to be 3a 2 /26 < 1. 

Figure 14.11 and Figure 14.21 illustrate the behavior of the empirical mean as the 
system transitions from having two equilibria to having a single one, which is con- 
trolled by the value of ^-m(O). This is an instance of a bifurcation of equilibria. 
From Proposition 12.21 we know that when h is small, the existence condition of two 
equilibria, ^-m(O) > 1, can be approximated by the condition 3a 2 /28 < 1. In the 
simulations, we let h = 0.1 so the approximate condition 3cr 2 /26 < 1 can be applied. 
In Figure I47T1 we change a but fix the other parameters, and consider the three cases 
^m(O) < 1 (3er726> > 1), |m(0) « 1 (3a 2 /29 = 1) and ^m(O) > 1 (3<7 2 /29 < 1). 
In Figure E3~2l we change 8. We can see that even though the parameters varied in the 
numerical simulations are not the same, the bifurcation behavior is similar. 

Figure [4~3l shows the effect of increasing h on the system stability. By stability we 
mean resistance to the transition of the empirical mean of the system from one state 
to the other (because the model is symmetric). The parameter h is proportional to 
the height of the potential barrier of each agent. Thus we increase the overall system 
stability if we increase the component's stability. This observation is analogous to 
comments in [311 ESI [26] . It is clear that h influences system stability substantially. 

Figure 14.41 illustrates the effect of system size on its stability. Clearly a larger 
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N=100, h=0.02, o=1, 9=10 N=100, h=0.1, o=1, 9=10 N=100, h=0.5, o=1, 6=10 

dm(0)/d^=1.0017, 3c 2 /29=0.15 dm(0)/d^=1 .0086, 3o 2 /29=0.15 dm(0)/d^=1 .0441 , 3o 2 ^8=0-15 




t 6000 a0D ° 10000 ° 0000 4000 t 0000 0000 10000 ° 2000 4000 t 



Fig. 4.3. The effect of changing h. Increasing it stabilizes the system. 



N=40. h=0.1, 0=1, 9=10 N=80, h=0.1, o=1, 9=10 N=1 20, h=0.1 , 0=1 , 9=10 

dm(0)/d4=1 .0086, 3c 2 /29=0.15 dm{0)/d$=1 .0086, 3o 2 /29=0.15 dm(0)/d^=1 .0086, 3o 2 /26=0.15 




Fig. 4.4. Influence of the system size N . A larger system tends to have a more stable behavior. 

system is more stable. These stability phenomena will be quantified with the large 
deviations analysis of Section [5j 

4.2. Heterogeneous Model. For the heterogeneous model, 9 is replaced by 9j, 
and the discretization is 

1 N 

X n+i = X n _ hU ( X f)At + ctAW 7 / +1 +Sj(—J2 X k- x ") At > ( 4 -2) 

fc=i 

with the same parameter settings. The different values of 9j are controlled by the 
parameters 9; and pi. In the simulation, we take K = 3 and {Qi}^ =1 = ©a/, } 
for a system a low, medium and high rates of mean reversion to the empirical mean, 
that is, the systemic risk. We also take {pi}fLi = {pl,Pm, Ph} for the corresponding 
fractions. We use the normalized standard deviation of the distribution of 9j values 
in order to quantify diversity. We find that the heterogeneous model behaves like 
the homogeneous one when h, a and ./V change. But, diversity on the rates of mean 
reversion has significant impact on system stability. 

As in the homogeneous case, in Figure l4~5l we consider cases with a below, close 
to and above the critical value. The results are similar to the homogeneous case 
as expected. For a below the critical value we have two equlibria and for a above 
the critical value one equilibrium. The condition ^m(0) > 1 is still necessary and 

sufficient for the existence two equlibria. The condition ^ j ^ 1 (/9;/0;)(3ct 2 /20; — 1) < 1 
is also a good approximation to the exact one when h is small. 

The parameter h and the system size N are closely associated with system sta- 
bility. We note that in Figure 14.61 and Figure 14.71 when h or N are increased, the 
system becomes visibly more stable. Another observation is that with h, a and N 
fixed, and with the mean of 9j of (|4.2|) equal to 9 of (|4.1[) . the heterogeneous system 
is consistently more unstable than the corresponding homogeneous model (see Figure 
14.31 and Figure I4.4[) . Clearly diversity tends to destabilize the system. 



9 



N=100, h=0.1, o=1, 0=[5;1O;15] N=100, fl=0.1, 0=2,1218, 0=[5;1O;15] N=100, h=0.1, o=3, 0=[5;1O;15] 

p=[0.33;0.34;0.33], stdev(6)/mean(e)=0.40825 p=[0.33;0.34;0.33], stdev(9J/mean(e)=0.40825 p=[0.33;0.34;0.33], stdev(e)/mean(e)=0.40825 

dm(0)/d^=1 .0096, E.(py0 1 )(3o72©-1)=-O,O949 dm(0)/d^1.0002, I.(p./e.)(3o /20-1)=1.O4O8e-O17 dm(0)/d^0.98922. I.(p./0.)(3a /20-1)=O.1219 




t ' t t 



Fig. 4.5. Effect of changes in a. The system has two stable equilibria when a is below the 
critical value and has single one otherwise. For small h, ^2fL 1 (pi /@t)(3cr 2 /2Qi — 1) < 1 is the 
approximate criterion. 



N=100, h=0.02, a=1, ©=[5;10;15] N=100, h=0.1, c=1, 0=[5;1O;15] N=100, h=0.5, o=1, 0=[5;1O;15] 

p=[0.33;0.34;0.33], stdev(e)/mean(e)=0.40B25 p=[0.33;0.34;0.33], stdev(e)/mean(9)=0.40825 p=[0.33;0.34:0.33], stdev(9)/mean(e)=0.40825 

dm(0)/d^=1 .001 9, E.(p./0.)(3o 2 /20-1 )— 0.0949 dm(0)/d^1 .0096, E.(p./0.)(3o 2 /20-1 )=-0.0949 dm(0)/d^1 .0494, Z.(p./0.)(3a 2 /20 r 1 )— 0.0949 




t t t 



Fig. 4.6. Effect of changing h. Increasing it stabilizes the system. 

We also change the diversity of 9j by changing 0/ and pi. To compare with the 
homogeneous case, in Figure l4~8l and Figure l4~9l we change the standard deviation of 
9j while the mean of 8j is fixed. In this most interesting part of the simulations we 
see that when we increase the standard deviation of diversity values, the number of 
transitions is notably larger than that in the homogeneous case. 

5. Large Deviations. In the previous two sections we saw both analytically 
and numerically that for large N, the empirical density X^{t, dy) is close (weakly, in 
probability) to the solution of the Fokkcr-Planck equation (I2.ip . and so the mean x{t) 
in stays around the first order moment of the deterministic limit, /_ yu(t, y)dy. 
If the condition of existence of two equilibria is satisfied, then x(t) will remain close 
to either — £b or +£b for relatively long time intervals, depending in particular on 
the parameter h. However, as long as N < oo, as we have seen in the simulations 
the random forcing by the Brownian motions {wjft)}^^ will cause transitions with 
non-zero probability. A systemic transition is the event that x(t) is displaced from 
to =F£b within a finite time horizon. Thus, systemic transition means that a 
large number of agents transition in a finite time. In this paper, we are interested 
in computing the probability of such a systemic transition. Mathematically, given a 
finite time horizon [0, T] and the conditions for existence of two equilibria, we want 
to compute the probability 

P(x(0) = -^,x(T)=^ b ) (5.1) 

when N is large and as a function of the parameters (h, 0, a) in ([l.lj) . 

5.1. Large Deviations of Mean-fields. According to [10], we can calculate 
this probability asymptotically for large N using large deviations. To state the large 
deviations theory that we will use, we will review briefly some notation and terminol- 
ogy from [TP] . 
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N=80, h=0.1, 0=1, ©=[5;10;15] 
p=[0.33;0.34;0.33], stdev(6)/mean(e)=0.41 059 
dm(0)/d^=1 .0096, E.(p/0)(3c 12® -1 ) =-0.0949 



N=120, h=0.1, o=1, 0=[5;1O;15] 
p=[0.33;0.34;0.33], stdev(6)/mean(e)=0.40907 
dm(0)/d|=1 .0096, E(p/0)(3c /20-1)=-O.O949 



N=160, h=0.1, o=1, 0=[5;1O;15] 
p=[0.33;0.34;0.33], stdev(e)/mean(6)=0.40758 
dm(0)/d^=1 .0096, Mp B ii3o^20-1)=-O.D949 




Fig. 4.7. Effect of changing the system size N . Larger system have a more stable behavior. 



N=100. [1=0.1,0=1, 0=[1O;1O;1O] 
p=[0.33:0.34:0.33], stdev(9)/mean(e)=0 
dm(0)/d^=1 .0086. E (p/0)(3o /20.-1 )=-0.085 



N=100, h=0.1, 0=1, 0=[6;1O;14] 
p=[0.33;0.34;0.33], stdev(9)/mean(6)=0.3266 
dm(0)/d^=1 .0092, Z (p./B )(3a /20.- 1 1 o 



N=100, h=0.1, 0=1, 6=[2;10;18] 
p=[0.33;0.34;0.33], stdev(6Vmean(e)=0. 
dm(0)/d^=1 .0089, E (p/0)(3a /20-1 )=-0. 




Fig. 4.8. The effect of changes in 0;. The median of the diversity values is fixed but the low 
and high sensitivities are changed to adjust the level of diversity of 6j while pi and the mean of 6j 
are the same. Increasing diversity tends to destabilize the system. 



• Mi (R) is the space of probability measures on R with the Prohorov metric 
p, associated with weak convergence. 

• C([0,T], Afi(R)) is the space of continuous functions from [0,T] to Mi(R) 
with the metric sup 0<t<T p(4>i(t), 4>2(t)). 

• Mr(R) = {fi G Mi(R"),7 <p(y)K d y) < where ip e C 2 (R) is a nonnegative 



function with lim 



kl- 



ip{x) = oo. From [10], if U(y) = y 3 



■ y, we can choose 

<p(y) 

• M oc (R) = U R>0 M R (R) = Mi(R), / ip(y)fi(dy) < 00} endowed with the 
inductive topology: fi n — ¥ fi in M oc (R) if and only if fi n fi in Mi(M) and 
sup n J (p{y)nn(dy) < oo. 

• C([0,T],M oo (R)) is the space of continuous functions from [0,T] to M oc (R) 
endowed with the topology: </>„(•) </>(•) in C([0, T], M^R)) if and only if 
M-) -> </>(■) in C([0,T],Mi(R)) and supo^y sup„ / (p(y)4> n (t, dy) < 00. 

• Given v G M OQ (R), we let = C([0,T] } M M (E)) : 0(0) = endowed 
with the relative topology. 

To simplify the notation, we rewrite (|2.ip as u t — C u u + hAi*u, where 



<9y 



y- / yip{t,y)dy 



M* 



dy 



[U(y)4>] 



Theorem 5.1. (Dawson and Gartner, 1987) Given a finite horizon [0,T], v G 
M 00 (R) and A C E u , if X N (0) = ± Y,f=i S x s (o) ^ v in M^R) as N ^ oo, then the 
law of Xjv(i) = jj 3xj(t) satisfies the large deviation principle with the good rate 
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N=100, h=0.1, o=1, O=[5;10;15] N=100, h=0.1, o=1, 0=[5;1O;15] N=100, h=0.1, o=1, e=[5;10;15] 

p=[0.1;0.8;0.1], stdev(e)/mean(6)=0.22473 p=[0.3;0.4;0.3], stdev(6)/mean(e)=0.38925 p=[0.45;0.1 ;0.45], stdev(e)/mean(6)=0.47673 

dm(0)/dS=1 .0089, I (p./0.)(3a /20.-1 )=-0.088 dm(0)/d^1 .0095, I.(p./0.)(3a /20-1 )=- 0.094 dm(0)/d^1 .0099, 1.(p./0.)(3a 720.-1 )=-0.0985 




Fig. 4.9. The effect of changes in pi, with ©; and the mean of 9j fixed. Increasing diversity 
tends to destabilize the system. 



function Ih-' 

- inf I h {4>) < liminf — \ogV{X N £ A) 

4>ak n^oo N 

< Km sup logP(X7v € A) < - inf 4(0), 
where A and A are the interior and closure of A in £ u , respectively, and 



i r 

I h{4>) = 7r^2 SU P J h{^,f)dt, 

to Jo f-.UJ^O 



(5.2) 



/oo 
f(y)<t>{dy), 
-co 

if 4>{t) is absolutely continuous in t £ [0,T] and Ih{4>) = oo otherwise. 

Remark. Here for <fi £ £" and t £ [0,T], </>(i) is viewed as a real Schwartz 
distribution on E, and A4*<j) are differential operators in the distribution sense, 
and / in (|5.2|) is a real Schwartz test function. The definition of absolute continuity 
for the path of measures (</>(i)) te [ 0i T] is hi the sense of Definition 4.1 in [TP] , that is to 
say: for each compact set K C K there exists a neighborhood of the null function 
in the set of test functions with compact support in K and an absolutely continuous 
function H K from [0, T] to R such that \(<j>(t),f) - (<p(s),f}\ < \H K {t) - H K {s)\ for 
all s,t £ [0,T] and f £ Ur- Note that by Lemma 4.2 in [10], if 4>(t) is absolutely 
continuous in t £ [0,T], <pt(t) exists in the distribution sense almost everywhere on 
t£[0,T]. 



In order to use Theorem 1 5. 11 wc let v = ul^ in (|2.2p and define the rare event A 
of systemic transition by 

A={^ef:^(T) = u|J. (5.3) 

However, since A is an empty set, Theorem 15.11 give a trivial lower bound for the 
probability in question. Therefore we consider instead the closed rare event A$: 

A s = {</> £ £ v : p(<f>(T),ul) < 5} . 

Then Theorem 15. II implies that 



inf I h {cj>) < liminf ^-\ogV{X N £ A s ) 

< limsup — logPpTjv £ A s ) < - inf I h ((f>). 
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In addition, we show that inf^g^ Ih{<P) can be bounded from below by ini^^ Ih(<P) 
as S — >• 0. 

Lemma 5.2. By definition inf^g^ Ih{4>) is decreasing with 5 > and bounded 
from above by inf^g^ 7/j (</>)• In addition, 

lim inf 7^(0) > inf Ih(4>)- 
Proof. See Appendix [Cl □ 

Combining Lemma 15.21 and the fact that inf , g ^ j 7/i(</>) < inf^g^ Ih{4>)t f° r an Y 
e > 0, we have for sufficiently small 5 > 

- inf J h (0) < lim inf i logP(Xjv € A«) 

< lim sup — logP(Xjv G A4) < - inf 4(0) + e. 

Therefore for large N and sufficiently small 8, 



P(X N £A 5 )nexp (-Nm£I h ((l>) . (5.4) 

This tells us that a larger system has a more stable empirical mean trajectory, which 
is consistent with what we have seen in the numerical simulation. Now the main step 
is finding inf^g^ Ih(4>), which is a min-max variational problem 

1 " T 



inf I h {4>) = inf — / sup {fa - C*^ - hM*<f>, f) 2 /(</>, f 2 )dt, (5.5) 
<peA 4>eA 2a l J / : (^,/a)^o v y 

where the / in the sup is a real Schwartz test function. 

5.2. An Alternative Expression for the Rate Function. The representa- 
tion of the rate function (15.21) is somewhat complicated, but we can simplify it if tj> 
has the density with some additional properties. If <j> is a density function such that 
</>(t,y) is smooth, rapidly decreasing in y 6 R for each t £ [0, T] and is absolutely 
continuous in t £ [0,T] for each y 6 R, then let g(t,y) satisfy 

<j> t -£%4>-hM*<j> = (<j>g) v . (5.6) 



Note that because of the properties of (/>, the left hand side of (|5.6p is well-defined in 
y G R and almost everywhere in t G [0, 7 1 ]. In addition, because <f> is positive valued, 
g exists and is unique except on a measure zero set in [0, T]. 
Note that for the pair ((f), g) satisfying 



sup J h W),f)= sup m),f y g) 2 /(<Kt),fi) = (<b(t),g 2 ), 

f--m)Ji)^o /:<0(t),/ y 2 >#o 

and therefore we have the following proposition. 

Proposition 5.3. If <p is a density function such that <f>(t) is a Schwartz function 
for each t € [0, T] and is absolutely continuous in t € [0, T] for each i/gl, and g(t, y) 
satisfies h5.b}) , the rate function Ihisf) * n (5.ty) can be written in the form 



W) = ^| (<P,g 2 )dt. 



(5.7) 
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We interpret (I5.6P and (|5.7p as follows. The function g is regarded as the driving 
force making <\> deviate from the solution of the Fokker-Planck equation (|2.1I) . and 
Ih{4>) is the L 2 ((j>) norm of g, which measures how difficult it is to have this deviation 

</>■ 

6. Small h Analysis. The goal of this section is to analyze the min-max problem 
(|5.5p which controls the asymptotic systemic transition probability. This problem is 
nonlinear and infinitely dimensional and is difficult to analyze. To get some useful 
information about it we will assume that h is small and analyze it in this regime. 
We will first solve (15 ,5|) when h is exactly 0, and then we will get rigorous upper and 
lower bounds for (|5.5p when h is nonzero but small. We will then compare the large 
deviations result with the local fluctuation theory of a single agent so as to explain 
why interconnectedness destabilizes the system. 

6.1. The h = and the Small h Analysis. We note that when h = 0, 

u ±& = u ±?o> where 



20 

In this case, (|5.5p is solvable and the optimal path is a Gaussian, starting from u e _^ Q 
and ending in u e + ^ Q . 

Theorem 6.1. Let h = and define 

P e (t, P ) = ^exp(- (y -°; a (t))2 |, a'(t) = ^t-&. (6.2) 

Then p e G A is the unique minimizer for \5. 5|) and 

inf :7 Q (0)=/ O (p e ) = % 



Proof. See Appendix [TXT] □ 

We show next that (|5.5p is continuous at /i = 0. 

Theorem 6.2. There exists j(h) such that j(h) 



as h -> 



inf ^(^ 

0eA 



2£ 2 



a T 



<l{h)- 



3.3) 



W^e recall here that 



6 - Co + fc& + 0(^ 2 ), 6 



ff2_6_ 

20^ 



1 - 2(a 2 /2d) 
l-3(<7 2 /2#)' 



(6.4) 



Proof. See Appendix [TJJ and UJJ □ 

As it is stated we could replace £& by £ in Theorem 16.21 since £5 = £ + o(l) as 
ft — > 0. We will see in the next section (in Proposition [7|5|) that 7(ft) = 0(h 2 ). In fact 
we show this rigorously for the upper bound but only formally for the lower bound. 
Since — £,o + h£i + 0(h 2 ) we see that the term 2£ 2 /(a 2 T) contains the leading-order 
term and the first-order correction in the /i-expansion of inf^g^ Ih(4>)- 
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6.2. Large Deviations for the First Exit Time. In this subsection, we con- 
sider the rare event B of systemic transition at some time before T: 

B 5 = {<f> € £" : 3* € (0,T],p(4,(t),ul) < 5}. 

In other words, Bg = Ut£(p,T]Ag(t), where 

A s (t)={0e£» : P m),ul)<6}. 

We let B := B . We then have that 

Lemma 6.3. By definition inf^ e B 5 Ih(<l>) is decreasing with S > and bounded 
from above by inf^gB Ih(<l>)- In addition, 

lim inf I h (<f>) = inf I h (<fi) = inf h(<p), 
5->o0es 4 #eu fe(0l 3.]i(t) </>eB 

wftere A(t) := A (i). 

Proof. See Appendix ID.4I □ 

From Theorem 16.21 we see that in the sense of large deviations the probability of 
system failure at some time before time T is essentially the same as the probability 
of system failure at time T. 

Corollary 6.4. For any t\ < t^, there exists a sufficiently small h such that 
inf^ e A(ti) h{4>) > inf^ eA ( t2 ) h{4>)- Consequently, inf h (<t>) ~ inf^eA(T) h{4>) for 
small h. 

6.3. Comparison with the Fluctuation Theory of a Single Agent. To get 

a better understanding of the large deviations results we need to carry out a standard 
fluctuation theory for a single agent. We assume that Xj(0) — —1 for all j and that 
the Xj(t)'s are in the vicinity of —1 so that we can linearize (ll.lj) : 

1 N 

Xj(t) = -1 + Zj(t), x(t) = -l + z(t), z(t) = xJ2 z ^- 
For V(y) = \y A —\y 2 , zj(t) and z(t) satisfy the linear stochastic differential equations 



a 



N 
\ " 



dzj = — (8 + 2h)zjdt + Ozdt + adwj, dz — —2hzdt + — dwj, 

3 = 1 

with Zj(0) = z(0) = 0. The processes Zj(t) and z(t) are Gaussian and the mean and 
variance functions are easily calculated. We are especially interested in their behavior 
for large N. 

Lemma 6.5. For all t > 0, Ezj(t) = Ez(t) = and Varz(i) = ^(1 - e" 4 ' 1 *). In 
addition, Varzj(t) — > 2 (e+2h) ^ — e -2(0+2/t)t) as N ^ oo, uniformly in t > 0. 

From Lemma 1631 we see that a 2 /N and a 2 /2(9 + 2h) should be sufficiently small 
so that linearization is consistent with the results it produces. 

6.4. Increased Probability of Large Deviations for Increased 9 and Its 
Systemic Risk Interpretation. We have now the analytical results with which 
we may conclude that individual risk diversification may increase the systemic risk. 
Assume that a 2 /N and a 2 /2(6 + 2h) are sufficiently small and N is large. From 
Lemma 16.51 the risk Xj(t) of the agent j is approximately a Gaussian process with the 
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stationary distribution J\f(— 1, <r 2 /2(9 + 2h)). If the external risk, a is high, then in 
order to keep the risk xj (t) at an acceptable level, the agent may increase the intrinsic 
stability, h, or share the risk with other agents, that is, increase 9. Increasing h is in 
general more costly (cuts into profits) than increasing 9, and at the individual agent 
level there is no difference in risk assessment between increasing h and increasing 9. 
Therefore the agents are likely to increase 9 and reduce individual risk by diversifying 
it. Note that a 2 /2(9 + 2h) < a 2 /29 when a 2 and 9 are significantly larger than h. 
Thus, individual agents can maintain low locally assessed risk by diversification, even 
in a very uncertain environment. 

What is not perceived by the individual agents, however, is that risk diversification 
may increases the systemic risk while it reduces their individual risk. Because a 2 and 
9 are significantly larger than h, the small h analysis can be applied and from (|5.40 
and Theorem 16.21 the systemic risk (the probability of the system failure) is 




for small 8 and h 
~2{a 2 /29) 



l-i{a 2 /29) 



0(h 2 ). 



We see that there are additional systemic-level a 2 terms in the exponent and £&, which 
can not be observed by the agents, increasing the systemic risk, even if the individual 
risk a 2 /29 is fixed. In other words, the individual agents may believe that they are 
able to withstand larger external fluctuations as long as their risk can be diversified, 
but a higher a tends to destabilize the system. 

7. A Reduced Large Deviations Principle for Small h. In Section I6TT1 we 
show that the large deviation problem inf^g^ Ih{(j>) is continuous in h so that we have 
the upper and lower bounds for inf^g^ Ih(4>) when h is small. In this section, we 
analyze with a formal expansion the optimal path for mf ( f > EAlh(4>) by assuming that 
it is of the form p e + O(h), motivated by the fact that the optimal path is p e for h = 0. 
In this way, we can obtain a reduced large deviations principle (a reduced Freidlin- 
Wentzell theory) for the systemic risk. That is, we obtain a reduced rate function 
corresponding to a finite dimensional system after ignoring higher order terms. The 
reduced rate function has all relevant information up to 0{h 2 ) terms, and we also 
need to expand cj) to 0(h 2 ). 

We assume that the optimal cj) — p + hq^ + h 2 q^ + . . ., where 

p(t,y) = -^L= exp \- {y = f ))2 1 , a{t) = {4>,y). 

In other words, we let the first moment of <fi be determined by a(t), and from the zero 
h case we know that ait) = a e {t) + 0(h). From the form of p and (|5.6[) , a natural 
parameterization for and q^ is the Hermite expansion 

n=2 y n=2 y 

Note that by the properties of p and a(i), (q^ , y n ) = (q^ , y n ) = for n = 0, 1 so we 
can start the Hermite expansion from n = 2. 
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The formal expansion result of this section is that if the optimal <j) — p + hq^ + 
h 2 q( 2 \ then 

-T / j 2 \ 2 

inf/ fc (0)w inf — / (-a + h(a 3 + 3^-a-a)) dt, (7.1) 

0SA a(t):0<t<T 2<7 2 7o \* 2 # ) 

a(0)=-£ 6 

for small /i. Note that a(i) = (0,2/) = x(t). The right hand side of (17. ip is an one- 
dimensional variational problem that has the form of a rate function of the Freidlin- 
Wentzell theory. In fact, the right side of (|7.1|) is the large deviations variational 
problem for the rate function of the small-noise stochastic differential equation 



dx(t) 



dt + eadw(t) (7.2) 



where here e = 1/\N is small. Note that 3cr 2 /29 < 1, as assumed above, and therefore 
(17. 2[) also represents a bi-stable structure. In the remainder of this section we describe 
how this result is obtained by formal expansions and then in Section [7.3l we show how 
we recover from (|7.1|) the main result of the paper stated in the previous section. 

An important remark about the expansion is that the Hermite functions are a 
basis of the L 2 space and thus p+hq^ +h 2 q^ is generally a signed measure. However, 
if g*- 1 -* and q^ can be expressed as the linear combinations of finite Hermite functions, 
then we can see that for any e > 0, there exists a sufficiently small h such that the 
negative part of p + hq^ + h 2 q^ is less than e. 

7.1. Optimization over g. The first step in finding the optimal <fi = p + hq^ + 
fi 2 q( 2 ) is determining the optimal g by using (|5.6j) for <fi. Once we obtain g, we can 
compute Ih{4>) by using (|5.7|) . It is also natural to assume that g = g^ + hg^ + h 2 g^ 
along with the Hermite expansion: 

.^^-E^iU 9 W = P -^Mt)^-P, ^=P- 1 E^)^-P. 

n=0 * n=0 " n=0 " 

In addition, since (cp~> ,y) = (q^ 2 \y) = 0, we can see that <p = p + hq^ 1 ' + h 2 q^> 
satisfies 

C%4> = C* pP + hC* p q {1) + h 2 C v q {2 \ M*(b = M*p + hM*q {1) + h 2 M*q {2) . 
The force U(y) = y 3 — y can also be expanded in Hermite polynomials: 

U{y)=p- l jZ5 n {t)^-p. 

n=0 y 

Now everything is expanded in the orthogonal basis and we can find the optimal g^ 
and g^ 1 ' by putting everything into (I5.6|) and comparing coefficients. 

Lemma 7.1. With the expansions mentioned above, the optimal g(°~) is —-^a, and 
the optimal (3 n for g^ are 

{p,U(y)}, n = 0, 

n+l)b n+1 -S n , l<n<3, (7.3) 
9(n+l)b n+1 , n>4. 
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Proof. See Appendix [EH □ 

It remains to determine g( 2 \ From (|5.7[) we see that the only contribution of g^ 
to Ih up to 0(h 2 ) is (p, 2g^g^} — — 27 ^a. Thus it suffices to determine 70, which 
can also be obtained from (|5.6[) . 

Lemma 7.2. With the expansions mentioned above, the optimal 70 is 

7o = -(g (1) ,t/(y)+g (1) ). 



Proof. See Appendix IE.21 □ 

7.2. Optimization over 0. We are now ready to find the optimal (j>. For given 
4>=p+ hqM + h 2 q^ and the corresponding optimal g = # (0) + + ft. 2 .g (2) , (|5J| 
gives 

' T ((p,(9 {1) ) 2 +2g^gW) + (q^,2g^g^)) dt + 0(h 3 ). 



h_ 

2^2 



From Lemma [7T21 (p,2#(°V 2 )) = -2g^(q^\ U{y) + .g (1) >, and therefore 

< Pl 2 ff (°y 2 )) + <g« 2<,«V 1) > = -2g^( q ^\U(y)) = -2g^£H n 8 n b n , 

n=2 

where H n (t) := (p"\ (d n p/dy n ) 2 ). We note that 

OO 

<p, 2 3 (°) 3 ( 1 )) = -2^*0, <p, 0,W) 2 > = 51 + £ <p, ( 5 ^) 2 ) = (, 9 (°)) 2 . 

71=1 

Then Ih{4>) can be written as 

4(0) = A / T (5 (0) - h8 Q fdt + £ f \h 1 (3 2 - 2H 2 g^5 2 b 2 )dt (7.4) 



+ £2/ (H 2 (3 2 -2H 3 g^S 3 b 3 )dt+^J2 [ H n p 2 n dt + 0^). 
°~ Jo a n=3 Jo 

We see that a and b n are coupled at the 0(h 2 ) level of (|7.4|) . However, from the 
results of the zero h case, a = a e + O(h) and p = p e + O(h) so we can decouple a and 
b n and express the expanded Ih(4>) U P to Olh 2 ) as the sum of independent terms. 

Proposition 7.3. To order 0{h 2 ), the rate function Ih{4>) can be written as the 
sum of independent terms: 

= A A.9 (0) - hS fdt + £ [ T (Hj 2 + 2^-a e H 2 5 2 b 2 )dt (7.5) 



2a 2 J ^ u; 2a 2 J y ^ dt 

h 2 r T d h 2 °° r T 

— J {H 2 ~P 2 + 2-a e H 3 5 3 b 3 )dt + _ £ jf HJldt + 0{h% 



n— 3 * 
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where H n (t) = ((p e )-\ (O n p e /dy n ) 2 ) , U(y) = (p 6 )" 1 E„=o |jrf> e > and 

-~6 = -(p e ,U(y)), n = 0, 

f t b n+1 + 9(n + l)b n+1 - ~6 n , l<n<3, (7.6) 
f t b n+1 + 6(n + l)b n+1 , n>4. 



We can see from (|7.5p that q^ does not appear in terms up to 0(h 2 ). From 
the h expansion of u±^ b in (|2.2j) . and the fact that V(y) is a polynomial of degree 
four, we have b n +i(0) = b n +i(T) = for n > 4. The variational problem for b n +i 
is to minimize J H n /3 2 dt where /3„ is given in terms of b n +\ by (|7.6p . The obvious 
solution of this problem is 6 n +i = and /3 n = for n > 4. Consequently, in order to 
find the optimal <f> for Ih(4>) in (|7.5p . we may solve separately the variational problems 
for a, bi, 62 and 63. 

7.3. Probability of Systemic Transitions for Small h. We consider the 
small probability of systemic transitions for large N and small h through the large 
deviation inf^g^ Ih(4>). Here we consider the solution up to O(h) terms. That is, 
using (|7.5[) . we solve the variational problem for a(t): 

inf / (gW-h5 ) 2 dt= inf / (^-a + h(a 3 + 3^a - a)) 2 dt. (7.7) 

a(t):0<t<Tj Q ^ ' a(t):0<t<Tj V dt V 26 " V ' 

a(0) = -£ b o(0)=-6 
a(T)=£ b a(T)=£ 6 

By simple calculus of variations methods we find the optimal a. 

Lemma 7.4. The optimal a(t) for J 7. 7| ) satisfies the second order ordinary dif- 
ferential equation 

^a = h^ + (3^-l)a)(3a 2 + (3^-l)) 
wii/i a(0) = — £{, and a(T) = £5. Consequently, the optimal path is 

a(t) = ^t-£ b + 0(h 2 ). (7.8) 

By inserting ()7.8[) into (|7.7|) we obtain inf <p e A Ih{$) U P to 0(h). 

Proposition 7.5. For small h, the large deviations problem, inf^g^ Ih(4>)> up to 
0(h), is 

inf 7,(0) = + 2^i) + 0(h 2 ), (7.9) 

where = £0 + + 0(h 2 ) from 112. 5\) . Note that £1 is positive because 28 > 3a 2 . 
Proof. See Appendix lET3l □ 

The asymptotic probability of systemic transition for large N and sufficiently 
small 5 and h has the form 

P(X N e A 5 ) « exp f-JV inf = exp f-JV {^(£0 + 2fc£i) + 0(^ 2 )]) . 



19 



8. Effect of Diversity of Sensitivities on the Transition Probability. We 

consider the situation introduced in Section [3] and analyze it when h = 0. We aim 
at computing the transition probability in this situation. The K partial empirical 
averages 



x k (t) := k=l,...,K (8.1) 



then satisfy a closed system of stochastic differential equations 

dx k — —Z=dw k {t) - 9 k (x k - x)dt (8.2) 
VPkN 

where w k are independent Brownian motions and the empirical mean x(t) can be 
expressed in terms of the partial averages as 

K 

x(t) = y^pfcXfcft) 
fc=i 

Proposition 8.1. If x k (0) = — £ b for all k = 1, . . . ,K, then x(T) is a Gaussian 
random variable with mean — and variance a\ := Var(x(T)) given by 

°$ = Jf£ Q T e Ms R- 1 (e Ms ) T Qds (8.3) 

where g is the K -dimensional column vector (p k )k=i,...,K , M and R are the K x K 
matrices defined by 

Mij = —9i(Sij — pj), R l] =p l 5 lj , i,j = 1, . . . ,K, 

and T stands for the transpose. 
Proof. See Appendix lF.il □ 

We can then deduce that the transition probability is 



Pt ~ exp 

> '/rp 



(-5) < M » 



Our next goal is to study the impact of the diversity on the transition probability. 
Proposition 8.2. Let us assume that the diversity is small: 

9 k = 8(l + 5a k ), <5«1 

where ^ fc p k a k = so that 9 is the mean value of the 9 k 's. The equilibrium position 
£b, the variance a\ and the transition probability px can be expanded as powers of S 
as 

k 



On 



N 



1 + f (l-e- S n 2 ds)+0(5 3 ) 

k •'° 



k 
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Proof. See Appendix lR2l □ 

This proposition shows that the diversity reduces the gap between the two equi- 
librium states and enhances the fluctuations of the empirical mean. Both effects 
contribute to the increase of the systemic transition probability. 

9. Summary and Conclusions. The aim of this paper is to introduce and 
analyze a mathematical model for the evolution of risk in a system of interacting 
agents where cooperation between them can reduce their individual risk of failure but 
increase the systemic or overall risk. The model we use is a system of bistable diffusion 
processes that interact through their empirical mean, a mean field model. We take the 
rate of mean reversion to the empirical mean 8 as a measure of cooperation, the depth 
of the bistable potential ft as a measure of intrinsic stability of each agent, and the 
strength of the external random perturbations a as the level of uncertainty in which 
the agents function. Using the theory of large deviations we calculate the probability 
that the system will transition from one of the two bistable states to the other during 
a time interval of length T, when the number of agents N is large and when h is small. 
In this regime of parameters we find that systemic risk increases with cooperation. 
The formula from which we draw this conclusion is given is Section 16.41 We also 
show that when the rate of mean reversion to the empirical mean varies among the 
different agents, that is, when there is diversity in the cooperative behavior then the 
probability of transitions increases, which means that the systemic risk increases. 

Acknowledgement. This work is partly supported by the Department of En- 
ergy [National Nuclear Security Administration] under Award Number NA28614, and 
partly by AFOSR grant FA9550-1 1-1-0266. 

Appendix A. Proof of Proposition [2T2l 

For small h, we view u| as a perturbed Gaussian density function. Let p^(y) be the 
Gaussian density function with mean £ and variance a 2 /29, Y be the Gaussian random 
variable with the density pg, and rj = 2/er 2 . By using the expansion exp(— hrjV) = 
1 - hryV + h 2 ri 2 V 2 /2 + 0(h 3 ), we have 




Then we calculate m(£) as follows: 




£ + hrj{£EV(Y) - E[YV(Y)}} + h 2 V 2 {--^EV 2 (Y) + £{EV(Y)) 2 




h V ^EV y (Y) + h 2 rj 2 ^{E[V(Y)V y (Y)] - EV(Y)EV y (Y)} + 0(h 3 ) 

2 2 

h v ^-EV y (Y) + h?r??-Cov{V y {Y),V{Y)) + 0(h 3 ). 
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The compatibility condition = m(^) gives 

Ey y (y) - hr]Cov(V y (Y), V(Y)) + 0(h 2 ) = 0. (A.l) 
Assuming that = £o + ^£1 + 0(h 2 ), the 0(1) terms in (IA.1|) give 



2 2 
£o+3^o-£o = £o(£ 2 + 3^-l) 



Then £ = 0, ±1/1 - 3cr 2 /20 if 3cr 2 < 20, or otherwise £ = 0. In order to obtain the 
nontrivial result, we suppose that 3cr 2 < 20 and £0 takes ±yl — 3a 2 /29 in the later 
calculations. Note that EV„(F) = £ 3 + (3cr 2 /26» - 1)£ = 2/i£g& + 0(/i 2 ), and 



Cov(i-„(n. nr.) = E[i/(y)^(y)] + o(fe) = e[(V - Vmv ; - n] + 

rl„» 3,., I, 



= E[-r 7 - -r 5 + -r 3 i + 0(h). 
L 4 4 2 

Along with the identity £g + 3<r 2 /20 = 1, we have 

EV' 3 = £„ + 0(h), EY S = ( 1 + 4^ - 6 Q J fo + 0(A), 

Then Cov(F 3/ (y),F(F)) = 6(ct 2 /26») 2 (1-2ct 2 /26»)Co + 0(/i). The 0{h) terms in (IA~T1) 
imply a = 3r,(a 2 /29) 2 (l - 2a 2 /20)/^ O - 

Appendix B. Proofs in Section [3l 

B.l. Proof of Theorem 13.11 The proof contains three steps. 

B.l.l. Existence and Uniqueness of the Weak Solution of the McKean- 
Vlasov Equation. The existence and uniqueness of a probability measure valued 
process (u%(t), . . . , Uk (i)) that is a weak solution of the McKean-Vlasov equation 
(|3T2|) is guaranteed by [IE1 Theorem 2.11]. 

B.l. 2. Weak Compactness of the Empirical Process. By Prohorov's the- 
orem, it suffices to prove that the sequence {(Xjf, . . . ,X^)}^ =1 is weakly compact 
by showing that 

sup sup sup E[{X^(t,dy), \y\}] < 00, 

N l<k<K0<t<T 

which can be done by using the calculations similar to (Bl) and (B2) in [5]. 

B.1.3. Identification of the Limit. For a test function / e <5(R), we define 
= (f(y),X l N (t,y)) = E ieIl f{ Xj {t))l\ll\- By Ito's formula, 



dX$ = — — [-hU(xj)dt + adwj + &i(x - Xj)dt]f y (xj) + -a 2 f yy (xj)dt 
1 iez, 

K 2 
= (-hUf y + 9 ; ((y,^P/Aj v ) - y)f y + ^-f yy ,X l N )dt + (f y , ^- £ S Xj d Wj ). 

1=1 1/1 jeii 
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Then by the integration by parts, we write 

K 2 

dx l N = {(hux l N ) y - [et({ y , E pi*n) - y)x l N ] y + ^(X l N ) yy }dt - E {8 Xj ) y d Wj . 



z=i 1 11 jeii 



For simplicity, we prove the case that K = 2 and the general case is similar. We let 

rl,xn v^-> 
L N x ^-N 



X]^ xn x X 2 ^ xn denote the product measure on 



X^ xn x X 2 ^ xn (y u y 2n ) = Xj.it, yi ) ■ ■ ■ X^t, y n )X 2 N (t ,y n +i) ■ ■■Xpf(t,y 2n )- 



For a test function / e <S(R 2rl ), we have 

d(f,X' N xn x X 2 N xn ) = d{.f,X l N xn x X 2 N x yV +d(f,X 1 N xn x X 2 N xn }( 2 \ 

where (1) and (2) denote the first and the second order terms of d(f,Xjj xn x X^ xn ), 
respectively: 

n 

d{f,X 1 N xn x X%* n )M = 'E&dXh&yj) x X^ n -^ j x X 2 N xn ) 

3=1 

2n 



1 " 

d(/,^ x "x^ x ")( 2 > = t £ (/,dXi(t l!/i )x^( i ,ft)x4 x ('- 2 »^xXf n ) 
+ 5 E (/,^(i,%)x^(t, yfe )x^ xn x^ x( "- 2) ^) 

j : k—n-\-l 
1 n 2n 



J=l fc=n+l 



Note that for j ^ k, dX l N (t, yj ) x dX l N (t,y k ) = ^TJ2iez,( S ^(yj))j( S xi(Vk))kdt = 
-^{5{y k - y^Xhfayj^jkdt, and dX^it^yj) x dX 2 N {t,y k ) = 0. If we analogously 
represent the generator G^ x i,xn X 2,x n ^f of (f,xjj xn x X^j xn ) as 



then G^ 1)Xn ,, x » / -)• as iV -)• oo and G^ liX „ 2>Xn / = G, x» ^.J, the gener- 

ator of (/, u xn x where (ui, u 2 ) satisfying ()3.2|) . Then the limit of (-Xjy, X^) is a 

solution of the martingale problem associated to (|3.2[) . In addition, by [18, Corollary 
2.10], the solution is unique and therefore (X x , X^) (u\, U2) weakly as TV — > 00. 
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B.2. Proof of Proposition 13.21 All we need to show is that for small h, 
Jp-m(O) > 1 if and only if a < a^ iv , where m(£) is defined by (|3.4I) . We obtain 
by calculate j| / yu e ls .(y)dy. Note that = (2&i/a 2 )(J yuf^dy - O z i,i and 



26; 
^•2 



l/u? i£ <fi/ _ 1 ] + | H-J 



Zi & \ y u it d y - € 



>2 

On the other hand, we can also compute j^Zi^ by directly taking the twice deriva- 
tives of Zi^: 



20i 



26/ 



By comparing (|B.1[) and (|B.2I) . 

d f e , 26, 

^1 y y^ d y = -r 



2/ "ffdj/ - ( / yufedy) 



(B.2) 



Note that J yuf dy = 0, so ^m(O) = X);=i /°z(20//er 2 ) / y 2 uf dy. By using the same 
trick in the proof of Proposition 12.21 let pi (y) be the Gaussian density function with 
mean and variance a 2 /26;, K; be the Gaussian random variable with the density 
pi, and rj = 2/a 2 . Then for small h, Z^ = 1 + /i^EV^) + 0(h 2 ), and 



y 2 uifidy 



y 2 (l-hr)V + 0{h 2 )) Pl {y)dy 



Z^(EY 2 - fc»,E[jf v(y,)] + 0(^ 2 )) 

Elf + fc»j(Elf EV(Yj) - E[lfV(Y/)]) + <3(h 2 ). 



Therefore |m(0) > 1 if and only if p/(26 / /cr 2 )(Ey ; 2 E^(Y ; ) - E[lfF(Yi)]) > 0. 



Note that EY, 2 = a 2 /2Q h EV(Yi) = (3/4)(E}f ) 2 - (l/2)Ey, 2 , and E\Y^V(Yi)] = 
(15/4)(EY; 2 ) 3 — (3/2) (EF; 2 ) 2 . Then the sufficient and necessary condition becomes 



K / 2 \ 



1=1 



2 

26; 



B.3. Proof of Proposition [373l It is equivalent to show that X>fc=i Pl/®l — 
~Yld = i pi&i Ylj—i Pi I® 2 - First note that by the Cauchy-Schwarz inequality, 



vJ=l 



i=l 



z=i 



Then it suffices to show that 1 < X)/=i Pi®i Tld=i Pi/@i- Again by the Cauchy- 
Schwarz inequality, 



A' 



A 



A 



Pi 



A 



1=1 



1=1 



6 



6, 



i=i 



24 



Appendix C. Proof of Lemma 15.21 

It suffices to show the case that S = 1/n. For each n, let <fr n G A 1 / n , such that 
inf^ £j 4 1/n Ih{4>) < h{<t>n) < ^4>eA 1/n h(4>) + 1/"-; {h(4>n)} are bounded from above 
by inf^g^ Ih{4>) + 1 < oo. Because Ih is a good rate function, and by Proposition B.13 
of [18], compactness is equivalent to sequentially compactness in C([0, T], M OQ (R)), 
{</)„} has a convergent subsequence whose limit fa is in A. As Ih is lower 

semicontinuous, then 

lim inf 4(0) = ]iml h ((p n ) = liminf h(fai k ) > h(fa) > jnf ' 4(<W- 
Appendix D. Proofs in Section [6j 

D.l. Proof of Theorem 16.11 We prove it in three steps. The first step is to 
show that there exists a uniform lower bound of Io(fa, for all <fi G A. 
Lemma D.l. Ifh = 0, then inf 0eA / o (</») > 2Q/(o- 2 T). 
Proof. For any e A, a(t) denotes J yfa[t,dy). We observe that 

J h {fa) = sup (fa - C^fa f) 2 /(fa f 2 ) f > (fa - C%fa y) 2 , 

because (fa 1} = 1. Note that (fa,y) = 4r(fay) — 4xCL{t), and 

(£%fa y) = (\a 2 fa y + 9^- [(y - a{t))fa , y) = -6((y - a(t))fa 1) = 0. 
Then after taking the infimum over <p € A, we have 

inf Ufa) > inf -L ( T ( ^-a] dt = inf -L f ( £ a \ dt = ?-§-. 

(p&A ~ (jj^-A 2<7 Jq \dt J a(ty.o<t<T 2a 2 J \dt J a 2 T 

a(0)=-£ 
a(T)=£o 

The last equality is obtained by a simple calculus of variation with the optimal path 
a(t) = 2( t/T - ft. □ 

The second step is to show that io(p e ) = 2£ 2 /(ct 2 T). Then inf^g^ Io(fa) — 
2^q /(tr 2 T) and therefore p e is a minimizer for (|5 . 5|) . 

Lemma D.2. Ifh = 0, and 

a t \ 1 f (y-a e (i)) 2 l 2^o 

p e ft,y)^^=exp^- Vtf Wy , a e (i) = ^-ft, 

^Trfi I % J 

tfien p e G 4 and 7 (p e ) = 2£ 2 /(cr 2 T). 

Proof. By reading (|5.6I) with = p e and ft. = 0, we have = C*„p e + (p e g) y - One 
can easily check that C* c p e — and pf = — Pl^o- e (t)- Then we have g = — j^a e (t) 
and by (joTTf . 



W) = 




□ 

Finally we prove that for h = 0, the minimizer p e is unique. 
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Lemma D.3. For h = 0, p e is the unique minimizer for 15. 5\) . 

Proof. From the previous lemmas, we find that if is a minimizer then a(t) = 
f y<p(t, dy) must be a e (t), and / = — -^a e (t)y is a global maximizer of Jq(4>, ■)■ Then 
for any test function /, 4: Jq(4>, ~^i ae (t)y + e /) = at e = 0. By a simple calculus of 
variations, <j> satisfies the linear parabolic PDE: 

1 d d 

k = + 6gj {(y ~ a e (t))0] ~ ^« e (*)^> 

with the initial condition 0(0) = w!le D , and that implies the uniqueness of the mini- 
mizer, which is p e . □ 

D.2. Proof of Theorem 16.21 (Upper Bounds). Define the test function: 

We recall that from (|2.3p and (|2.5p . £b depends on h and — > £o as ft 0. 
Proposition D.4. For any e > 0, tften /or aZZ sufficiently small h, 

inf 7,(0) < -L /" V, (l a « ^ fc(„3 _ y)) 2 >ttt + £ (D 1} 

It is not difficult to see that the first term of the right hand side of (|D.ip is equal to 
/ ( a2 T) up to a term of order h as ft — > 0. 
Proof. We construct the test function 0" £ A as follows: 

f (1 - + *g[0,5T], 

<T(i) =^ u (t), te(5T,T-5T), 
l(l^)f"(i) + ^ b , t G [T — <5T, T], 

where oT will be determined later. Note that inf0 gj 4 Ih{4>) < Ifi(4> u ) so we just need to 
compute 4(0"). Let satisfy dES]) for = 0". For t G (<5T, T - 5T), (f> u {t) = p u (t), 
and it is easy to see that pj 1 = -f t a u p u y and £;„p" = 0. Therefore for f G {ST.T-ST), 
g u = -j- t a u - h(y 3 - y) by l[5T5]) . From (|5T7|) . we have 

The rest is to show that for any e > 0, there exists a sufficiently small h such that 
the last term in the last equation is bounded by e. It suffices to show that for any 
ST > 0, we can choose a sufficiently small h such that (<fi u , (g M ) 2 ) is bounded by a ST- 
independent constant c" > for t G [0, ST] U [T - £T, T]. If so, then let oT < ea 2 /c u 
and 

i(^ + £J <0 "• (! '" )2><i^s i (2<T)c " <e • 
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for sufficiently small h. 

For t G [0,<5T], because <j) u is simply the convex combination of uie 6 and p u , 4> u 
can be bounded by a <5T-independent constant. To compute g u from (|5.6|) . it is also 
easy to see that C^ u (f> u and A4*0" can be bounded by 5T-independent constants. The 
only term we need to worry is {p u (t) — u e _^ b )/5T from computing <pf(t). However, p u (t) 
is differentiable at t = and p"(0) —> u e _^ h as h — > so we can bound {p u (t) — u e _^ b )/8T 
by a 5T-independent constant with suitable /i. Thus is bounded independently of 
<5T and so we can find a ST- independent constant c" > such that (</>", (g u ) 2 ) < c u . 

The same argument works for t e [T — ST, T] and we have the desired result. □ 

D.3. Proof of Theorem [672] (Lower Bounds). From (|D.1I) . there exists some 
constant C such that inf^g^ Ih(4>) < C for all h < ho- Then we can assume that 
Ih{<t>) < C for all e A and all h < ho without loss of generality. The following 
lemma shows that the first and second moments of all <f> £ A are uniformly bounded. 

Lemma D.5. Given C > 0, there exists R > such that for any <f> G A with 
Ih(4>) < C f or some h>0, then 



Proof. Recall that Mr(R) = {<j> G Mi(R),/ tp(y)(f>(dy) < R} and Afoo(R) = 
U.R>oA/fl(R) with the inductive topology. Here we focus on the case that i/3 = 1 + y 2 
in order to obtain the uniform result, and let M|(R) and M^(R) denote the spaces 
with the quadratic Lyapunov function ip. 

The proof is an application of Theorem 5.1(c), Theorem 5.3 and Lemma 5.5 of 
[TD]. By Theorem 5.1(c), if G C([0, T], M^,(R)) with 0(0) = ul Cb and ^(0) < C 
for some /i > 0, then is in an ^.-dependent compact set K. By Theorem 5.3 the 
compact set K is contained in C([0, T], Mn(l)) for an /i-dependent i? > 0. Finally, 
by Lemma 5.5 and Theorem 5.1(c), it suffices to let i? > e AT (C + r), where r and A 
satisfy 



r>2 if(y)u e _ ib (y)dy, A> sup (fi, C^ip + hM<p + -<fil)/(v, f), 

with y>(j/) = 1 + y . Obviously we can find the uniform r and A for all h > and 
also the uniform R. Then any of interest are in C([0,T], Mj|(R)) and thus have the 
uniform bounded first and second order moments. □ 

Now we derive that lower bound. The key idea is that because we have the 
universal upper bound for the first and second moments of all <f> G A and for all 
h < ho, Chebyshev's inequality implies the uniform convergence. 

Proposition D.6. For any e > 0, then for all sufficiently small h, 



Proof. Define / = i * f , where / is a piecewise linear function and i is the 
standard mollificr: 



sup (<Hi),y) 2 < sup m),y 2 ) < R. 



te[o,T] te[o,T] 



1 




(D.2) 




otherwise. 
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Then f M is a smooth function with the compact support [— 2M — 1,2M + 1]. In 
addition, / M (y) = y on (-M + 1, M - 1), < 1, and \f™\ is uniformly bounded 
for all M and is nonzero only on U|__ 2 (iM — 1, zAf + 1). 

Because for all <f) G A, (0(t), (fy 1 ) 2 ) < 1, we can estimate the rate function: 

4(0) > ^2 / V " £ ^ - ^*<^ f M )* dt 

> ^ {j\<k - c;^ - hm*^ f M )dt 

Then we estimate the integrand term by term. By Lemma ID.5| the following conver- 
gences are all uniform in <fr G A and h < ho. 
First we have 

[ T (4, t ,f M )dt=(ul b J M )-(u^ b J M ). 
Jo 

u Mt, are ex P onen tially decaying functions so (u±^ b ,f M ) converges to ±£b rapidly as 
M -» oo. 

We note that (£^, f M ) = a 2 (0, /#)/2 - 6(4,, (y - a)/f ). By reading the prop- 
erties of fyy and Chebyshev's inequality, we have (0, f^) — > as M — > oo. We write 
(4>,(y-a)f™) as 

(0, (y - a)/ H M ) = a(l - (0, /*')) + ((0, y/f ) - a). 

Since a is bounded and (0, fy 1 ) -> 1 as M -> oo, a(l - (0, -> as M -> oo. We 
see that 

|(0,y/, M )-a| 2 < (2 / |i/|^(iy) > ) 

Y »/(-M+l,Af+l) c / 

< 4 / y 2 0(dy) / 

J(-Af+l,M+l) c J (-M+1,M+1)<= 

Again by Chebyshev's inequality, the right hand side vanishes as M — > 00. 
Finally we estimate (A4*4>, f M ). Since f M is compactly supported, 

\(M*4>J M )\ = 1(0, (y 3 - y)0\ < (2M + l) 3 + (2M + 1). 

For a fixed M, we can choose a sufficiently small ft, such that h\(J\A*<f>, / )| is small. 

Consequently, for any e > 0, we can first choose a sufficiently large M and then 
there exists a sufficiently small h such that 

mf : W) > % ~ e. 




□ 
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D. 4. Proof of Lemma 16.31 It suffices to show the case that S = 1/n. For each 
n, let 0„ e B 1/n and t n G (0,T] such that p((f> n (t n ) , u|J < 5 and Laf^ eBl/n 4(0) < 
4(0n) < inf0GB 1/ „ 4(0) + 1/"-; {4(0n)} are bounded from above by inf^gs 4(0) + 
1 < co. Let {t nk } be a convergent subsequence of {£„}. Because 4 is a good rate 
function, and by Proposition B.13 of [18], compactness is equivalent to sequentially 
compactness in C([0, T], Moo(R)), {0 nfc } has a convergent subsequence {0™ fc , } whose 
limit 0* is in where t* = lmit nk . As 4 is lower semicontinuous, then 

lim inf 4(0)=liminf4(0rv)>4(0*)> inf 4(0) > inf 4(0). 

n </>S-Bi/„ n fc ' 0eA(t*) 4>£B 

Appendix E. Proofs in Section 

E. l. Proof of Lemma 17.11 We note that p t — —p y -^a and therefore 

d , ^ d , d n , d ^ , d n+1 

n=2 y n=2 y 

, 2 ^ d Q n ,,d ^, <9™ +1 

+ * L i M Cn W P ~ dt^^dy^ 1 " 

n=2 a n=2 a 

After collecting 0(1) terms in (j5.6[) and integrating over y, we have 
-P^a= ^ 2 Py + ^(2/-a)p + TO (0) =W (0) - 

Then = -^o. 

Now we collect 0(/i) terms in (15.6[) and integrating over y. We get 

gn+l 



d d n d d n 1 

n=l n=2 y n=2 a 



P 



n=2 y 71=0 y 71=2 y n=0 J 



Using the fact that 



1 2 d n+1 „, ,9™ „ 



we have 

E ^ & *+^ = ~* 2> + + E S ^P + E A^*> 

71=1 a 71=1 " 71=0 3 71 = W 

and the optimal /3 n are obtained by comparing the coefficients. 

E.2. Proof of Lemma 17.21 Let V'' 2 ' denote the anti-derivative of that 
vanishes at — oo. After collecting 0(h 2 ) terms in (|5.6[) and integrating over y. We 
have 

^ (2) = iaV 2 > + % - a)qW + u{y)q^ + g< 2 V 0) + q^gW + pg™ . (E.l) 
2 J 

Note that pg^ — J2^=o Jn-§-^P, so 70 is obtained by integrating (|E.1[) from y = —00 
to y = 00. Then we have 70 = — ((r 1 ', t^(y) + <7^)- 
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E.3. Proof of Proposition We write a(t) = a (t) + ha x {i) + 0(h 2 ) with 
a (t) = 2£ t/T-£ and ai(t) = 2£ 1 i/T-£ 1 . Then we put a(t) into (|7!7|l and we have 

We note that is a constant, and ao(t) and £lo(t) are °dd functions with respect 
to t = T/2. Then 

1 f T f ( d \ 2 d d 



2io (e + 2/ia) + o(/i 2 ). 



Appendix F. Proofs in Section [8J 

F.l. Proof of Proposition!!]} The system of SDEs flgrj for the vector 
(5fc(i))fc=i,...,_ff has the form 

dX(t) = MX(t) + -^=R-^ 2 dW{t) 



where W(t) = (u>k{t))k=i k is a column vector. This system can be solved: 



X(t) = e Mt X(0) + -£= / e M ^Rr^ 2 dW(a) 
vN Jo 



If Xfe(O) = — £b, then, using the fact that the uniform vector is in the null space of M, 
we have e Mt X(0) = X(0). As a corollary we get the explicit representation of the 
empirical mean: 

x(t) =_& + -£= /' g^e^-^R-^dWis) 
v N Jo 

This shows the desired result. 



F.2. Proof of Proposition [8721 The expansion of £ 2 follows from the explicit 
expression ()3.5p . The expansion of o\ follows from the expansion of (I8.3|) and uses 
the properties of the matrix M. We have M = —9M — 58N, with 

M = I — ug T , where u = (1, . . . , 1) is the if-dimensional column vector, 
Nij = ai(Sij - pj), i,j = l,...,K. 

The matrix M satisfies M™ = M for all n > 1 and therefore 

e -6Mt = y t^M n = I + Y =1+ (e~ St - l)M. 

We have 



71! 1 — ' 71! 

n—0 n—1 



oc 



t- ^ 77 I 



n 
n=0 
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Using the fact that M T g = (and again that M™ = M for n > 1), we can expand 
g T e Mt ^ g T + Sg T{(-0t)N + (e- St - 1 + 6t)NM) 

+ 5 2 g T \^N 2 + [e~ 8t -l + 0t- @$-] [N 2 M - 3(NM) 2 + NMN] 
- 6t[e~ St -1 + et] {NMfg} + 0(6 3 ). 

Using the fact that M T N T g = N T g and M T (N T ) 2 g = (N T ) 2 g, this can be simplified 
into 

g T e Mt = g T + Sg T {e -9t _ 1)N + s 2 g ^[(et) 2 - (i + et)(e~ Bt - 1 + 0t)] N 2 + 0(5 3 ). 

Consequently 

g T e Mt R-\e Mt ) T g = g T (I + (e~ St - l)M)R-\l + {e-' et - l)M T )g 
+ 28g T (e- St - l)NR~\l + (e~ St - l)M T )g 

+ 25 2 g T [{et) 2 - (1 + M){ e - St - 1 + h)] N 2 R-\I + (e- St - l)M T )g 
+ 5 2 g T {e- St - l)NR-\e- et - 1)N T g + 0(S 3 ). 

Using the fact that M T g = and NR~ 1 g — Nu — 0, we obtain 

g T e Mt R-\e Mt ) T g = g T R- 1 g + S 2 (l - e^f g T N R- 1 N T g + 0(S 3 ) 

We have g T R~ x g — 1 and g T NR~ 1 N T g — J^k Pk a 1 which gives the expansion of the 
variance a\. 

Finally the expansion of the transition probability can be obtained by substituting 
the expansions of Q and a\ into (|8.4[1 . 
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